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Abstract 

According to one-loop perturbation theory, fermions whose masses 
are totally generated from Yukawa couplings do not decouple in the 
heavy mass limit. We investigate this issue nonperturbatively in a 4- 
dimensional Zi scalar-fermion model with staggered fermions. Our data 
at intermediate and stronger Yukawa couplings on 8'^ and 12^ lattices 
suggest the nondecoupling of heavy fermions as predicted from one-loop 
calculation. However, at the strongest Yukawa coupling where a possible 
multi-critical point may come into play, we cannot be conclusive. 

1 Introduction 



The decoupling theorem says that when a particle has a mass much higher than 
the physical scale, it will have very small influence (e.g. radiative corrections) on 
the "physical world". However, it can be easily seen at one-loop level that the 
contribution of the particle to renormalized quantities of other particles will not be 
suppressed by its huge mass if its mass is generated from the Yukawa coupling, i.e. 
if its mass is generated through the mechanism of spontaneous symmetry breaking 
(SSB). This is the so-called phenomenon of nondecoupling. 

Nondecoupling of heavy fermions in theories with SSB has been discussed in 
several papers f^, H, ^. Until last year, all arguments were within (one-loop) per- 
turbation theory (except in [^] where large-N expansion was used). Nondecoupling 
beyond one-loop was still not clear and should be explored. 

Recently, we studied this issue in a nonperturbative way in a ?7(1) scalar-fermion 
model with explicit mirror- fermions There, we ignore the gauge field assuming 
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that this approximation will not change the picture qualitatively. Hence, the scalar 
field is our "physical world" . The relevant quantity to measure to decide decoupling 
or not is the ratio of the two renormalized Yukawa couplings: One is the usual 
renormalized coupling defined as the renormalized ( mirror- )fermion mass divided 
by the renormalized vacuum expectation value (VEV); the other is obtained from 
the fermion-fermion-scalar 3-point vertex function. In the weak coupling regime, 
one-loop calculation shows that the ratio of the two is very close to one. As the bare 
Yukawa coupling gets stronger, the renormalized Yukawa coupling defined from the 
mass-to- VEV ratio is no longer a good definition for the renormalized coupling. 
There, we need to define the renormalized coupling from the appropriate 3-point 
vertex function. If the ratio of the two Yukawa couplings (denoted by R) still 
remains 1.0, we take it as the indication for the nondecoupling of the heavy (mirror- 
)fermion. 

Our main results in are: 

(i) Up to the strong Yukawa coupling regime, the ratio R for the mirror-fermion is 
always equal to 1.0 within the error, indicating that the heavy mirror-fermion does 
have its renormalized Yukawa coupling proportional to the mass and will not decou- 
ple from the transverse component of the scalar field. Thus one-loop picture survives 
the strong Yukawa coupling limit. The idea of decoupling the mirror partners by 
giving them large masses does not work. 

(ii) Since the action of the model has a symmetry between fermion and mirror- 
fermion , it is obvious that the heavy fermion itself will not decouple either. 

(iii) Our data also show that all doublers decouple as expected. 

(iv) We think that these conclusions also apply to the SU{2) version of the mirror- 
fermion model because its qualitative behaviour appears to be similar to that of the 
f/(l) model If. 

(v) We conjecture that in other scalar-fermion models, heavy fermions do not de- 
couple either. 

In order to confirm our conjecture mentioned in (v), we carried out Monte 
Carlo simulations on a 4-dimensional Z2 scalar field theory coupled to the stag- 
gered fermion. We now report on our results in this letter. 



2 Action and Renormalized Quantities 

The lattice action of the 4-dimensional Z2 scalar-fermion model with staggered 
fermions is 

s = -2Kj2Y.<f>^ <^.+^ + A J2{<f>i - 1)2 + ^ 02 

X /1=1 X X 

1 ^ 

+ 7;Y.Y. Xx ?7x,M {Xx+^, - XxJ + 5'int (1) 

^ X fj.= l 

where (p^ is the scalar field with only one real component, Xx and Xx are the one- 
component staggered fermion fields, r/^,^ being the phase factor associated with 
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staggered fermions, rj^^i = 1 and r]x,ii = (— for 2 < /i < 4, 5'int represents 
the term for the Yukawa interactfon. In this report, we study the model with the 
so-called "non-overlapping" Yukawa coupling introduced in |^. Hence, 

Sint = GY^Xx^xXx , = 

X x' Ghypercube{x) 

where the set of sites x' is obtained by 

-X 



+ ^ti , = 0, 1 , 



and [.] denotes the largest integer part of the arguement. The boundary condition 
is chosen such that it is periodic for the scalar field and the spatial directions of the 
fermion fields, and antiperiodic along the temporal (4th) direction of the fermion 
fields. Since each species of staggered fermions will generate four Dirac fermions in 
the continuum with opposite chiralities, the physical spectrum will be vector-like. 
The above action has the following discrete symmetry: 

0x ^ -<Px , Xx ^ , Xx • (2) 

For the study of the decoupling of heavy fermions, the relevant phase (called the 
FM phase from now on) is the one where Z2 symmetry in eq.(^ is spontaneously 
broken. Notice that we do not have the bare fermion mass term in the action. This 
means that the renormalized fermion mass will be totally generated by the Yukawa 
coupling in the FM phase. Definitions of various renormalized quantities there are 
as follows. 

The renormalized scalar mass m/j and the scalar wavefunction renormalization 
constant are defined as 

Z 

^2 = „ . {4>{pS*{Pi))c (3) 

P4 -|- TTlj^ P4^P4min 

where subscript c means the connected part of the propagator, 
^(P^) = ^T. , Pi = 4sin2(|^ 
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Q is the total number of lattice sites, x means (x, X4), and 4'*{pa) is the complex 
conjugation of 0(^4). In our simulation on a finite lattice, we set p^mm = 27i/L. 
The renormalized VEV is defined as 



VR = ^ , 101 = 0, . (4) 

\JZ^ " X 

Similarly, the renormalized fermion mass fiR and the fermion wavefunction renor- 
malization constant Zp are defined as 

J2e-^^^y^{XyXo) ^ ZFiTnip,))-' (5) 



fiR + i sin(p4) 
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where the spatial 3- momentum is set to zero, is chosen to be tt/L for the fermion 
fields, and Tjilp^) is the renormalized fermion 2-point vertex function in momentum 
space. The renormalized Yukawa coupling Gr is defined as 

Gn^^. (6) 

The relevant quantity to tell decoupling from nondecoupling of heavy fermions 
in this Z2 scalar-fermion model of staggered fermions with nonoverlapping Yukawa 
coupling is still the ratio of the 3-point renormalized Yukawa coupling and the one 
defined in eq.(||). The 3-point renormahzed Yukawa coupling coupled to the scalar 
field is defined as 

Gg^ = -M±^ TM G(^) fn{q.) , (7) 



where k^, p^, are the 4th components of the momenta of scalar field, fermion and 
anti- fermion, respectively, and A;| is 4 sin^(k4/2). We have set the spatial components 
of all momenta to zero. The appearance of the Kronecker-delta above is due to 
energy-momentum conservation. G^^^ is the connected part of the (j)-x-X 3-point 
Green's function and is 

In our simulations on lattices we choose 

271 71 71 

The ratio R is defined to be i? = G^^ /Gr. 

Notice that SSB only occurs in the zero-momentum mode, the above connected 
3-point Green's function G^'^'> is equal to the disconnected one, because it is defined 
aX k = (0, ki) with k^ 7^ 0. In the Z2 scalar-fermion models, there are no massless 
Goldstone bosons in the FM phase since no continuous symmetry is spontaneously 
broken. We simply treat the scalar field itself as the "physical world" and assume 
that the issue of the decoupling of heavy fermions does not depend on whether 
the broken symmetry is discrete or continuous. Notice that refiection positivity for 
scalar-fermion models with staggered fermions cannot be proven. We assume that 
no ghost particles are present in the spectrum. 



3 Phase Structure and Numerical Simulation 

The phase structure of the model defined in eq.(0) at A = 0.01 has been explored 
and was presented in figure 2 of 0. At G = 0, the model reduces to a one- 
component pure scalar theory which has FM, symmetric (PM) phases and an 
anti-ferromagnetic (AFM) phase. At A = 0.01, the transition point between FM 
and PM phases occurs around k = 0.120 and is a Gaussian fixed point. Due to 
the symmetry: (p^ — > (— 1)^0^,k —Ik, the transition between PM and AFM 



4 



phases occurs around k = —0.120. As G is gradually turned on, the FM-PM phase 
transition occurs at smaller and smaller (and eventually negative) k values because 
dynamical fermions favour the FM phase. This feature has been observed in all 
scalar-fermion models studied so far 0, §. As the value of G is further increased, 
some new phase, like the ferrimagnetic (FI) phase, shows up and the phase structure 
will be (slightly) dependent on the model 0. For the Z2 scalar-fermion model 
investigated in this letter, the FI phase and a new PM phase were found at strong 
and infinite G Around G = 1.7, there is a point (called point B in [|^, and will 
be called as such from now on) around which four phases (FM, PM, AFM and FI) 
may coexist. We have explored the phase structure around that point B on 4^ ■ 8 
lattice, and find that point B is located around G = 1.7 and k, = —0.125, consistent 
with previous estimate 0, |^ . 

The physically interesting phase transition is the one between the FM and PM 
(called SI in 0) phases (curve AB in figure 2 of 0). It was found to be consistent 
with a second-order phase transition on which the cutoff can be removed. We 
therefore carry out Monte Carlo simulations at A = 0.01 in the FM phase, but in 
the vicinity of this FM-PM phase transition line. 

In the simulations, besides A = 0.01, we set G to be 0.3, 0.6, 1.0, 1.2 and 1.8. 
The lattices we use are 8^ and 12^. We then tune the scalar mass parameter n to 
have rrifi around 1.0 on 8^ lattice to reduce finite size effects and at the same time 
maintain a reasonable cutoff. 

In this Z2 scalar-fermion model, no continuous symmetry is spontaneously broken 
in the FM phase, so there is no massless particle in the spectrum. Finite size effects 
should be small except on small lattices where vacuum tunnelings happen. It is well 
known that on a finite lattice, the ground state is nondegenerate. This means that 
there is no SSB on a finite lattice. The unique ground has a vanishing VEV and is 
a symmetric linear combination of two states, which are peaked at the positive and 
negative minima respectively. It is these two states that will converge to the two 



degenerate ground states in the thermodynamic limit [T^. Once tunnelings happen 
on a finite lattice, we will be in the nondegenerate ground state of the system. 
Finite size effects will be dominated by an instanton-like equation [jlO| and may not 



be small on small lattices. If the system is not too close to criticality such that it 
is trapped in one of the minima, then finite volume effects will be dominated by 
perturbative effects rather than tunneling. In the presence of dynamical fermions, 
the precision of our data will not be good enough. To analyze finite size effects 
dominated by tunnelings is therefore too demanding for the moment. So, when we 
study the properties of the FM phase of the model on a finite lattice and eventually 
extrapolate to get the infinite volume limit, we actually do not wish to see tunneling 
events. 

Without loss of generality in our simulations, we always set up the initial condi- 
tion such that the system starts in the positive minimum. We define the time-slice 
average of the scalar field as 

X 

By observing values of those (psit), we will know whether tunnelings happen or not. 
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The Monte Carlo simulations are performed by the unbiased Hybrid Monte Carlo 



method ||TT|. Therefore, the fermions have to be doubled by taking the adjoint of the 
fermion matrix for the second species. (The fermionic part in eq.(|Tp is given for a 
single species of staggered fermions.) We will have eight degenerate Dirac fermions 
in the continuum limit. The number of leapfrog steps per molecular dynamics 
trajectory was chosen randomly between 3 and 10. The step size was tuned to 
maintain an acceptance rate around 75%. The necessary inversions of the fermion 
matrix were done by the conjugate gradient iteration, until the residuum was smaller 
than some small value times the length square of the input vector. We find that 
this value has to be 10^^^ on the 8^ and 12^ lattices. We use Creutz observable 
exp{~6H) to decide whether the system has equilibrated or not |12| where SH is 



the difference between the new and old Hamiltonians in the Hybrid Monte Carlo 
update. In equilibrium, we should have 



(e-^^) = 1 . 



4 Conclusion and Discussion 

Our numerical data are presented in table 1. Our data on the renormalized fermion 
mass and bare VEV at G = 1.2 and 1.8 agree with those published in 0. We also 
find that the renormalized fermion mass has the mean field behaviour //^ = 
The renormalized fermion mass presented here can be as high as 700 to 800 GeV at 
strong Yukawa couplings. (The physical scale is set by vr = 246 GeV.) At all points 
where we did simulations, no vacuum tunneling events were observed. We think this 
is partly due to the fact that dynamical fermions tend to increase the height of the 
barrier between the two minima. It is also because Hybrid Monte Carlo is a local 
updating algorithm. (In the presence of dynamical fermions, we cannot use cluster 
algorithms to perform global updates.) Thus it is not easy for the system to tunnel 
even on an 8^ spatial lattice. 

According to table 1, data on R from G = 0.1 to 1.2 are all consistent with 
1.0 within errors, indicating a universal behaviour of the nondecoupling of heavy 
fermions as expected. 

However, at G = 1.8, k = —0.04, data on R on 8^ lattice (i.e. point g in table 1) 
is clearly smaller than 1.0 by several standard deviations, while data on 12*^ lattice 
are too noisy to tell. Therefore, we cannot be certain whether heavy fermions are 
still coupled to the scalar field at G = 1.8. At the moment, we would like to say 
that there could be a possibility that heavy fermions do decouple as the decoupling 
theorem says, although their masses are totally generated from Yukawa coupling. (In 
the action of our model, we do not have a bare fermion mass term.) This possibility 
exists because points f, F, g, G in table 1 are lying slightly to the right of point B 
and may very well be affected by it. Since point B might be a multi-critical point 
and may have, in principle, renormalization properties different from those dictated 
by the Gaussian fixed point at G = 0, it should not be too surprising that heavy 
fermions decouple in the vicinity of point B. Besides, if we stay at G = 1.8 and keep 
reducing the value of k, we will be approaching the phase transition line between 
the FM and FI phases. This FM-FI transition line is the place where bare VEV is 
still nonzero. Thus, it should not be a physically relevant phase transition line. 
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At this point, we would like to conclude that at least along the phase transition 
line between the FM and PM phases from G = 0.1 to G = 1.2 where the system is 
still governed by the infrared stable Gaussian fixed point at G = 0, heavy fermions 
whose masses are totally generated by the Yukawa coupling do not decouple. One- 
loop picture is qualitatively correct throughout this region. Although data presented 
in this letter are at A = 0.01, we believe that the above conclusion holds at all values 
of A between zero and infinity. As we go to even stronger Yukawa coupling where a 
possible multi-critical point may come into play, the decoupling of heavy fermions 
remains a possibility. However, our present data cannot give a conclusive signal for 
decoupling. Whether heavy fermions really decouple or not there depends on the 
properties of that possible multi-critical point. Due to limited computer resources, 
this issue will be left for the future study. 
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Table 1: Our Monte Carlo data are presented here. Points with lower-case letters 
are obtained on 8^ lattice while points E and F are on 12^ lattice. Each point has 
about 30000 molecular dynamics trajectories for equilibration, and around 100000 to 
300000 for measurements. But point G has only 30000 trajectories for measurements. 





G 




Vr 


mn 


f^R 


Gr 


^(3) 
^R 


R 


a 


0.1 


0.129 


0.600(3) 


0.578(8) 


0.115(1) 


0.191(2) 


0.205(43) 


1.07(9) 


b 


0.3 


0.128 


0.942(3) 


0.928(11) 


0.525(2) 


0.557(2) 


0.546(60) 


0.98(7) 


c 


0.6 


0.110 


0.654(2) 


0.90(1) 


0.709(4) 


1.083(4) 


1.10(10) 


1.02(8) 


d 


1.0 


0.080 


0.496(2) 


1.27(3) 


0.876(2) 


1.766(6) 


1.88(17) 


1.06(9) 


c 


1.2 


0.0(30 


0.137(2) 


1.15(3) 


0.907(2) 


2.08(2) 


2.07(21) 


0.99(9) 


E 


1.2 


0.060 


0.430(2) 


1.34(14) 


0.906(3) 


2.11(2) 


2.17(66) 


1.03(20) 


f 


1.8 


0.00 


0.362(2) 


2.29(6) 


1.142(2) 


3.15(4) 


2.97(44) 


0.94(5) 


F 


1.8 


0.00 


0.368(4) 


2.21(54) 


1.140(2) 


3.10(8) 


2.2(1.2) 


0.71(40) 


g 


1.8 


-0.04 


0.296(1) 


1.37(3) 


0.766(2) 


2.59(2) 


2.15(24) 


0.83(5) 


G 


1.8 


-0.04 


0.275(7) 


1.47(36) 


0.768(7) 


2.79(10) 


2.9(1.4) 


1.03(46) 
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